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ABSTRACT 


This thesis investigates the mechanical behavior of a hard 
polyurethane rubber in compression in the temperature range of 0°F 
to 70°F. A compression test is chosen because data on this type of 
deformation are lacking and the material is often used in compression. 
The neo-Hookean, Mooney and Ogden models are shown to be inadequate 
in predicting the elastic response of this material. A new model 
based on the Rivlin-Saunders form of the strain energy function is 
proposed. An attempt is made to model the viscoelastic behavior of 
this material with a simple nonlinear Kelvin-Voigt model. The effect 
of temperature on the elastic behavior is insignificant in the 
temperature range studied but the viscoelastic response appears to 


be affected. 
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CHAPTER I 
INTRODUCTION 


The subject of this investigation is the determination of 
the mechanical behavior of a polyurethane rubber at various temperatures. 
The material is assumed to be homogeneous, isotropic, and incompress- 
ible and thus may be treated as a continuum. The theory of finite 
elasticity and viscoelasticity is used to develop a theoretical model 
for this material. 

Previous investigators of rubber-like elasticity have used 
tension tests and shear tests to provide data. Since there is a lack 
in data from compression tests and also because the material under 
investigation is often used in compression, a compression test is to 
be used. Use of such a test will eliminate the stability problem of 
necking encountered in tension tests. 

One of the problems involved in this investigation is the 
temperature control system. To determine if the material behavior 
is temperature dependent, the sample must be maintained at the test 
temperature within a tolerance of + 2°F. Another problem arises 
from the nature of the compression test; that is, the friction be- 
tween the platens of the subpress and the sample during a test, as 
shown in Figure 1(a). 

The theoretical problem is to derive a model that will pre- 


dict the response of the material. The material is assumed to be 
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FIGURE: 1 LOAD DIAGRAM AND COORDINATE SYSTEM 
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a nonlinear Kelvin-Voigt solid. If the samples are loaded as shown 
in Figure 1, it is assumed that a spring-dashpot combination, shown 
in Figure 2 may be used as the mechanical model. It is assumed that 
the spring represents the elastic part of the stress and the dashpot 
represents the viscous or dissipative component of the stress. | 

A rectangular Cartesian coordinate system shown in Figure 
1(b) is chosen to facilitate the use of Cartesian tensors. The sum- 
mation convention is followed unless otherwise noted. 

The particular material studied is Vibrathane 6006, a poly- 
ester-based urethane elastomer, which has been cast to a hardness 
of 55 to 56 on the Shore D scale. The specimens are cured for 3 
hours at 312°F and then 12 hours at 160-170°F. The oil resistant qua- 
lities of this material make it highly suitable for oil field service. 
In one application, the material is subjected to large compressive 


loads in temperature range of O°F to 80°F. 
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CHAPTER II 
THEORY 


Assumptions 


The material is assumed to be a homogeneous, isotropic, in- 
compressible solid which may be modelled by the nonlinear Kelvin- 
Voigt mechanical model. The basic element of this model is a spring 
in parallel with a dashpot. With this assumption, the stress tensor 
may be decomposed into two parts: an elastic part and a viscous part. 
It is assumed that the elastic part may be derived using the methods 
of finite elasticity; the viscous part is derived from the theory of 
viscoelasticity. In the derivations of theoretical models, the effect 


of temperature is neglected. 


The Elastic Component of Stress 

The constitutive equation for a purely elastic solid, repre- 
sented by a spring, may be derived mechanically without consideration 
of thermodynamic effects. The stress components are defined as single- 


valued functions of the deformation gradients: 


g@) = fF) (1) 
where 5 (e) = the elastic stress tensor 
f = a tensor function 
F = the deformation gradient tensor 
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Since rigid body rotations do not affect the stress, the function, f, 


must satisfy 
Q £(F) Q' = #(QF) (2) 
where Q is any proper orthogonal tensor. From polar decomposition 
Peaeku (3) 


uniquely; R is a rotation tensor and U is a positive-definite, symmetric 


deformation tensor. Thus equation (2) may be written 


£(F) = Q! #(QRU)Q (4) 
Now since Rr! is an orthogonal tensor, substitution of Rr! for Q in 
equation (4) gives 
a il, 
f(F) = RF(U)R (5) 


~~ ~~ Dw 


which shows f to be a function of the six independent components of U. 
Also it may be shown that this is a necessary and sufficient condition 
for equation (2) to hold. Due to the difficulty in obtaining R and 


U, let 
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Then 


o(°) = Fh(y) FT (7) 
Introducing the Green deformation tensor, C = ue = FIF, 
o(C) 5 h(U) (8) 
and equation (5) becomes 
og) = Fg(c) FT (9) 


~ ~ 


From the assumption of isotropy, ¢ must be an isotropic function of 


b= ol + 440 + O90" (10) 


where dg: d1> bo are scalar functions of the scalar invariants of C and 


I is the unit tensor. The invariants are 


Loy EriGs 
Ip = 5 [ltr ¢)* - tr C4] (11) 
I, = det, C 
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Alternatively, the invariants and equation (10) may be expressed in 


terms of the Cauchy deformation tensor, B = Fae giving 
(e) _ 2 3 
ON TEE are tO? SE 
Using the Cayley-Hamilton theorem, 


Bo - 1,B¢ + 1 


B 18 geoplal = 0 


equation (12) becomes 


2 
ae = Vues, Voe 


where the response functions Vo» Vy> bo are new functions of I|; 


13. From the condition of incompressibility, 


equation (14) reduces to 


gl) = - pl + 4B + VB 


(12) 


(13) 


(14) 
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(15) 


(16) 


where p is a Lagrangian multiplier representing a arbitrary hydrostatic 


pressure and ~., w> are functions of I,, I, only. Thus equation (16) 
eae Pe 


is the constitutive equation for a homogeneous, isotropic, incompress- 


ible elastic solid. 
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The use of thermodynamic concepts results in a more restric- 
tive form of equation (16). A strain energy function, which deter- 
mines the mechanical properties of the material, is postulated. This 
function is equivalent to the free energy for isothermal deformations 
and to the internal energy for isentropic deformations. Assuming that 
the strain energy is a single-valued function of the deformation 


gradients: 
W = W(F) ; Ci.) 


an energy balance will yield 
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(18) 
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Comparison of equation (18) with equation (16) shows that the response 


functions of equation (16) have the form 
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In Cartesian tensor notation equation (18) becomes 
(e) OW OW _ 2 Ww 
Sig. = TaPy Ge ep ome, eee WKCKS (20) 


For simple uniaxial compression as shown in Figure 1, it 
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is assumed that the sample is free to expand radially and that no 
barrelling occurs. Let the deformation in the axial direction be 2, 


given by 


Q 
A=> (21) 
2g 
where 2 is the deformed length of the cylindrical sample and Ly is the 
undeformed length. For the coordinate system shown in Figure 1(b), 
this corresponds to the principal extension ratio, h- From the in- 


compressibility condition, 


hyAgdg = 1 (22) 
the other extension ratios are 
hy = d, = 4 (23) 
In terms of the extension ratio, the strain invariants are 
LPs he See ee 
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The Cauchy deformation tensor is 


ONO 
Ba 1g Heli (25) 
Ee IN 

0 0 4 


Since Con = 03 = 0 and Oi5 = 0 for i # j, the only component of 
stress remaining is 071° Letting O11 = Bele the constitutive equation 


(20) for this deformation reduces to 
ae (26) 


The form of W is to be determined empirically. 


Forms of the Strain Energy Function 

Before equation (26) can be used a suitable form of the 
strain energy function must be found. In determining the form of the 
function, the mechanical properties of the material are formulated. 
Material symmetries impose restrictions on the form of W which, in 
general, is a function of the three principal extension ratios, 
Ay» hos and Ag. Isotropy restricts the form of W to symmetrical 
functions of the principal extension ratios. Thus W may be expressed 
as a function of the three strain invariants given by equation (24). 
The condition of incompressibility restricts W to a function of I, 


and I, only. 


A majority of the previous theories in finite elasticity 
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assumes that the most general form of W for an incompressible, iso- 
tropic, elastic material is 
5h EO yh DOs ekg at ee w (27) 
iz0,j-0 J! : 
with Cog = 0 at zero strain. The simplest form of equation (27) is 


the strain energy function for the neo-Hookean material, 
I,-3) (28) 


For small strains the constant Ci reduces to E/6, where E is Young's 
modulus, and Hooke's law is recovered. The Mooney material has a 


strain energy function which is dependent on I,: 
W= C1 (14-3) + Cy; (15-3) (29) 


Similarly, for small strains, Cio + Coy = E/6 giving Hooke's law. 

These two simple forms of W are proposed as approximations 
to the strain energy function of the material under investigation. 
More complex forms use additional terms from equation (27) but the 
complexity of the resulting constitutive equation is increased. A 
common problem with the proposed strain energy functions is the lack 
of generality of the constitutive equations. Many strain energy 
functions yield constitutive equations that give good correlation 
with experimental data for a limited range of deformation or a 


limited number of types of deformation. 
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Of the more complex forms of W proposed in the past, one 
of the more promising and simpler versions is that suggested by Rivlin 


and Saunders [11]: 


W = C, (1-3) + (1-3) (30) 


where the function f is to be determined empirically. This form is 
obviously a generalization of equations (28) and (29) but more re- 
strictive than equation (27). Another form is that proposed by 
Ogden [9]. Finding equation (27) too restrictive, Ogden begins 


with the more general form: 
W = W(Aq sAp 0g) (31) 
and then forms the new strain invariants, 
J(a) = (Ay MHA Ag) 3) /0 (32) 


where a iS a real number. These invariants are related to I, and I, 


by 


Cc. 
— 
ine) 
~~ 
iT) 


] 
De (1,-3) 
and J(-2) = - 5 (15-3) 


Ogden proposes a strain energy function which is a linear combination 
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of these new strain invariants: 


Wrestle (34) 


where C. are constants related to the classical shear modulus, wu, and 


Young's modulus, E, by 
Co,.=2.=> (35) 


The Viscous Component of Stress 
The viscous component of stress, represented by a dashpot, 
is a dissipative component which in general is a function of the de- 


formation and the rate of deformation: 
o!) =~ pI + g(B,d) (36) 


where d is the rate of deformation tensor. Expanding 


z 2 2 
g(B,d) = Kol + KB + K,B° + Kad + Kyd 


(B2a+dB") 


drag Bt) a take ( Bie 


+ Kk, (Bd“+d°B) + Ky(B d@+d°B*) (37) 


where K.(1=1,2,3,...8) are polynomials in the joint invariants of B 
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14 
and d. Assuming that the viscous stress is a function of d only, equation 


(36) becomes 
go) =~ pr + k(d) (38) 


with k(0) = 0. Again from the assumption of isotropy, k(d) is an 


isotropic tensor function of d and thus may be written 
k = yal + yad + Yod° (39) 
Soi) Otty dives 12° 


where Yo? Y7> and Yo are scalar functions of the invariants of d and 


hans 0 if the invariants of d are zero. Then 
gV) = - pl + Yol + tyd + Yo" (40) 


which is the constitutive equation for a Reiner-Rivlin fluid. 


Equation (40) reduces, for a linear viscous fluid to 


oY) = - pl + Etr d I + 2nd (41) 


~ 


where & and n are constants which are independent of d. If the fluid 


is also incompressible, tr d = 0, and thus 


g\V) = - pl + 2nd , (42) 


~ 


n is the viscosity coefficient. In Cartesian tensor notation, 
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equation (40) may be written: 


(v) _ : 
Outen Oak GgEPlO, tH ved 41 tayo dapdys (43) 


In the case of simple uniaxial compression 


ees). 
Siler an (a) 


The condition of incompressibility gives 


det =50ys (45) 

from which it follows that 
Pea able as A oe 
Since spel? = o446¥) = 0; doo = d34 = - \/2hs the remaining viscous 


stress component, 0°) = oid in an incompressible linear viscous 


fluid under uniaxial compression is 


o!¥) = 3n O (47) 


The viscosity coefficient, n, is to be determined empirically. 
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The Nonlinear Kelvin-Voigt Solid 

The basic element of the Kelvin-Voigt model consists of a 
spring in parallel with a dashpot. In the linear theory both compo- 
nents are linear. A nonlinear model results when one of the com- 
ponents is nonlinear. 


The response of the linear spring to a force F is 


F = kx (48) 


where k is the spring modulus and x is the extension. Similarly, 


the response of the linear dashpot is governed by 


where n is the viscosity. 

In this investigation a creep test is used to determine 
the viscoelastic response of the material. A load is applied, held 
constant, and the deformation is recorded as a function of time. The 
response of a linear Kelvin-Voigt model is governed by the dif- 


ferential equation 
kx + nx = CH(t) (50) 


where x, the deformation, is a function of time, C is a constant, and 


H(t) is the Heaviside unit step function. 
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From the initial condition x = 0 at t = 0, the solution of equation 


(50) is 
C k 
x(t) = aR - exp (- aah H(t) (51) 


The generalized Voigt model consists of Kelvin-Voigt elements in 
series with a Maxwell element, which is a spring in series with a 
dashpot. The response is 

n 


t) H(t) + C Py) 


ke 
{1 - exp(- ate t)} H(t) (52) 
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thy et 
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where Ka and Ny are the constants of the Maxwell element and n is the 

number of Kelvin-Voigt elements. The first term, representing the 

response of the Maxwell element, allows for instantaneous elasticity. 
For the nonlinear Kelvin-Voigt model it is assumed that 

the total stress is the sum of the elastic component of stress and 


the viscous component of stress: 


(e) + o(V) (53) 
Thus combining equations (14) and (40) gives 


= Vol + 4B + VB" Spiel tae Y9d° (54) 


The assumption of incompressibility reduces this to 


2 
g=- pl + ¥B + vB + yd + y7d oe 
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which is the constitutive equation for an isotropic incompressible 
Kelvin-Voigt solid. 

The simplest nonlinear Kelvin-Voigt model is one where the 
spring is nonlinear and the dashpot is linear. Then for the case of 


uniaxial compression, equations (26) and (47) give 
(56) 


Further reduction is obtained if a neo-Hookean material is assumed 


for the nonlinear spring: 


AAS Pree (57) 


o = 20, 0A x 


gh 
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Now if a step load, CH(t), is applied, equation (57) may be written 


CH(t) 2 


CHE) = ac, (ae - +) + an 4 (58) 


LE 
rv 
where A = Ao/* is the current cross-sectional area of the sample. 


Letting K = CH(t)/2A the variables can be separated in equation 


01° 
(58), yielding 
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In general, the cubic polynomial, Ng - Kx? - 1 may be separated into 
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(A-R) (A°+PA+Q) , where R is the real root of the polynomial. Equation 
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(59) is then integrable and the solution is a transcendental equation 


having the form 


2 
en 37 5 fin es : AERP) ta pel sa Zr + 5 (60) 
AC, (QRP4R™) (X+PA4Q)  ¥(4Q-P*) 7(4Q-P4) 


where S, the constant of integration, is 


. 2 

= 3n (1+R) 2 - 2+P 

Ss = ———; | In (61) 
4C, (Q-RPAR*) fn Csr 7aa-r%) "Jaq. rea 


This equation is used as a first approximation to the viscoelastic 
behavior of the polyurethane material. Inspection of equations (56) 
to (59) shows that numerical methods are required for more complex 


forms of W and if n = n(A). 


If the Mooney form of W is chosen for the spring, the con- 


stitutive equation for the elastic part becomes 
1 2 
- x) (C, + a) : (62) 


Numerical integration is required to solve the resulting differential 


equation. The form suggested by Ogden [9] yields 


ote) = Gy eget), (63) 
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and when combined with the viscous tensor also requires numerical inte- 


gration, unless r = 1 and d= 2 in which case the neo-Hookean form 
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is recovered. Similarly, the form suggested by Rivlin and Saunders 
[11] in its simplest variation gives the same result as the Mooney 


form and other variations give rise to complex nonlinear differential 


equations. 


20 


as BO 


tartnorstttb er xatanos. 07.8 sot ovig epeiteriol ot 


a 
i 
+ 
a4 
~ 
2 
3 
= of 
as See 
v 
oh 
m 
Sy 
= 
i 
ie 
se $32 ; = 
ae 4 
s , s As 
i 
a ( i 
S 2 
is r 
a 
if 
} 
7 4 
“ 
fe 
4 = 
i 
, _ 
e 
“y “ 7 - 
yy) 
» } V 
a I 
~~, 1 
> ‘ae 1 y 
- 
= é - — 4 el yy 
Lo « 
i 
us = j 
. eA 
*. =—— 7 
a ss +e x 


CHAPTER III 
EXPERIMENTAL INVESTIGATION 


Description of Apparatus 

The apparatus consists of two separate systems: the en- 
vironmental control system and the compression testing system. Basi- 
cally, the enviromental control system is a refrigeration system which 
is capable of lowering a sample from room temperature (70°F) to 0°F 
and maintaining a set temperature within a tolerance of + 2°F. The 
compression testing system consists of a subpress and a testing 
machine. A schematic drawing of the apparatus is shown in Figure 3. 

The refrigeration system, except for the evaporator coil, 
is supplied by a commercial refrigeration manufacturing firm. It 
consists of a single stage 1/3 horsepower reciprocating compressor, a 
water-cooled condenser, a drier, a solenoid valve, a thermostatic 
expansion valve, and a thermostat control. The compressor is air 
cooled by a fan. A solenoid valve is placed upstream of the ex- 
pansion valve to minimize "overshooting" of the desired temperature. 
This valve closes when the compressor shuts off and prevents re- 
frigerant from passing through the evaporator. A 2-1/2°F temperature 
differential thermostat is used to control the compressor. The 
sensor is placed on the subpress block to take advantage of its 
thermal inertia and the fact that its temperature variation approxi- 


mates that of the test specimen in the platens. The refrigerant is 


ra 


i) — 


Pea) pa 
\ yorTantTzayat, JNTVENTAISAXS ar hg 


i 
{ se 
J 4 7 
i - ' 2 
E x ii ore ey re 
UIE = 4 oR Ay tae 
; 

; : 


“Ne SHI semodeyve adarsqpe ‘Ows: to edaredo eud'pied 


-teek 9322 pitt2as dragons sey bns mageye Fordiioo. laa" | 
oF neat donerentaten: a, med 22 gorabs) isdeonnbtvne at's 

40 ob G08) stud eradiias inoor moyt atamse 6 prtiene > ol 

ont 4°S + to sons-9fod 5 Wis Fw pberegne? Jaa brfintd : 

outteod 6 bas ores ‘5 to 2t2ztanos mnoeyerontie not! ie 

£ sword nf mwode 2h -ausstsqq6 std to prtwarb: at ino? a & nti | 

efi Foo noteronsys att 10% dijsone? ematere nottenaptvten at ia a ) 

ry mvt? ont hud asttism. Apt ievoptyiay ‘eign aNd: strane. at 

“6 NogzeavgHOD pai tenorgtoen NenafoeTot £\h sapeta’ oferta. i 40° 

si dst 2omvany. ix ovtsv bhonalde 6° 9Pb 6: J atigbaos viheecsll a 

; its at wos2antines) oat <feaetiag tateomany 6 bie ov lev ae a 

a eee baie to mesma2qs, boost ‘si aviay’ Bronoice A * tet Nd bofons 

; _srtenactos Absit zob- ont to" "ot sontetever asiainia of avisv noteneg * a 

a atnavarg bas. To. vedi siete sAs nant eseots SVT ov att “ : 

| smut evoqnss ide S$ A Wot erogsvs att fipuorAs, pnizzeq wich famed i 

a ooh ia ~ozasranes ants foystiag 03. mts at Yedzom ‘edd Tate tgm 


= 


i 


ay 


re A) 


- =" Ee .” - 2 > Pe 
oc - sad : xe Belpre. ms se) . i a - - 
a 8 <a 7 ha L oe ry ‘ if 7 a] 


KEELY 


MOOnnwnrd- ODAN OA FWON — 


| 


Compressor 
Condensor 

Drier | 

Solenoid Valve 
Expansion Valve 
Subpress Block 
Compression Cylinder 
Evaporator Coil 
Light 

Upper Platen 
Sample 

Lower Platen 
Thermostat Sensor 
Metal Inserts 


Access Window iF plate > baal @ 
‘ rR 


(a) Refrigeration Unit 


Lower Compression 
Head 3 : 
Upper CLD AMM PATI ITA VELMA : 
Compression : 
Head 
Polyethylene 
Shroud 
Thermostat 


S 
y 
@) 
‘ 


FA. 
VIR EM OE J QURSEN G) 


iPVELA 


Z) 


as 
LA 


(b) Environmental Chamber 
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freon 12. 

The evaporator coil of the refrigeration system is the cool- 
ing coil of the environmental chamber. This chamber consists of two 
concentric cylinders of galvanized iron with the cooling coil soldered 
to the inner cylinder in the wall space between the cylinder. Heat 
transfer from inside the chamber is thus improved and styrofoam in- 
sulation is used to reduce transfer from the outer shell. The base 
of the chamber is a 3/4 inch thick Socder plate with four 1/8 inch 
diameter steel inserts embedded in it to support the subpress. The 
top, also a 3/4 inch thick wooden plate, has a circular opening in 
it to accommodate the compression cylinder of the subpress. A poly- 
ethylene shroud covers the exposed shaft of the compression cylinder 
and the opening to prevent the formation of ice on the shaft which 
renders it immovable at temperature below 32°F. Access to the inside 
of the chamber is made through a doorway in the mid-section of the 
cylinder and a plexiglass door is fitted to allow viewing. Inside, 
the subpress is attached by fins to two rings, which are in contact 
with the inner wall of the chamber, for greater heat transfer. With 
this arrangement, it is possible to maintain the subpress temperature, 
and thus the sample temperature, within + 2°F of the desired temper- 
ature. Below the access doorway is a shelf which provides storage 
space for the specimen. A small light is placed above the doorway 
for illumination. 

The subpress component of the testing apparatus is a Tinius 
Olsen model which has been modified for these particular testing con- 


ditions. Modifications include the revision of the lower platen, 
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which was free to swivel, to make it rigid and parallel to the top 


platen; the addition of fins to provide a more efficient heat transfer 
path. The testing machine used is made by Gilmore Industries. It 

has a 50,000 pound capacity hydraulic fatigue test frame with a 
modular electronic control panel. The compression piston has a maxi- 
mum travel of 6 inches. For the tests the lowest load range (0 - 5000 
pounds) is used. The load to the nearest ten pounds and the displace- 
ment to the nearest thousandth of an inch are read directly from a 
digital voltmeter in the control panel. 

Miscellaneous equipment include a potentiometer, a switching 
box, a strip chart recorder, and a micrometer. Copper constantan 
thermocouples with an ice bath as the reference junction are used to 
measure the temperature. A thermocouple is embedded in a sample which 
is placed between the platens of the subpress. Then under test con- 
ditions, a position is found on the subpress block where the temperature 
and its variation closely resembled that of the test specimen. The 
sample with the embedded thermocouple is left on the shelf to indicate 
the temperature of samples sitting on the shelf. Another thermo- 
couple is attached to the subpress block at the position located by 
the trial and error procedure mentioned above to indicate the temper- 
ature of a sample in the subpress. These thermocouples are connected 
to a switching box which is connected to the potentiometer or the 
strip chart recorder. A micrometer, which is accurate to the nearest 
thousandth of an inch, is used to measure the dimensions of the samples. 


The design specifications of the environmental chamber and 


subpress modification are given in Sheriff [12]. 
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Preparation of Samples 


The polyurethane material presents some difficulties in 
the preparation of test samples. A commercial manufacturer provides 
the fetenia in 3/4 inch diameter rods which are 3-1/2 inches long. 
These rods are cut into 3/4 inch and 1 inch lengths on a lathe using 
a spring collett chuck. A coolant is used in the cutting operation 
to prevent excessive temperature rise in the sample. This gives a 
reasonable finish as well as parallel ends. To achieve as smooth 
a finish as possible, the ends are polished by successively finer 
grades of abrasive paper on a belt grinder. A Surfmet grinder, model 
16-1255, made by Buehler Ltd., fitted with Speed Wet Durite Cloth 
belts is used with water as a coolant. Final polishing is done by 


hand on very fine abrasive paper. 


Experimental Procedure 


The polished ends of a sample are lubricated with a silicone 
lubricant and the sample is placed between the compression platens 
of the subpress. The compression heads of the testing frame are 
applied to the subpress compression piston but with zero load. The 
chamber, subpress and sample are allowed to reach the desired test , 
temperature. Then the sample is removed and a micrometer is used to 
obtain initial dimensions of the sample at test conditions. Centering 
the sample in the subpress platens again, the sample is allowed to 
reach equilibrium Barer applying a load. The first load increment 
is applied and held until the deformation reaches a steady constant 


value, which is recorded with the load and the temperature. Also the 
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time required for the deformation to reach its equilibrium value is 
recorded. These steps are repeated for further increments of load. 
After the final load increment has been applied, the sample is un- 
loaded completely, measured, and then allowed to recover in the 
chamber at the test temperature. By making measurements of the sample 
length at regular intervals, a recovery curve may be obtained. The 
equilibrium load-deformation data give the stress-strain curve and 
thus shows the static elastic response of the material. 

To determine the viscous response a creep test is performed. 
A step load is applied and held on the specimen and the resulting de- 
formation is recorded at regular time intervals. This test gives a 
creep curve for the material at the test temperature and repetition 
of the test at other temperatures gives the effect of temperature on 


the viscous response. 


Analysis of Data 

The static elastic response data, obtained by allowing the 
deformation to reach equilibrium, are shown in Figure 4 in the form of 
stress versus extension ratio curves. Four different temperature 
tests are shown indicating the effect of temperature on the elastic 
response of the material. Also plotted on this figure are the various 
responses predicted by the theoretical models. The neo-Hookean and 
Mooney predictions are inadequate for X < 0.95. A single term ver- 
sion of Ogden's suggested strain energy function is also inadequate 
in predicting the material's static elastic response. 

Since the simple theoretical models are inadequate, a new 


model based on the Rivlin-Saunders form of the strain energy function 
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is proposed. It is assumed that the stress-strain curve is composed 
of two elementary curves arising from the two terms of the strain 
energy function. The first term of equation (30) by itself results 
in the neo-Hookean model. By an appropriate choice of C, the neo- 
Hookean model correlates with the experimental data for i < 0.5. 
Then to fit the data for 0.5 < > < 1.0, the second curve must have 
the form K(1-A"), where a > 1 and K is a constant. Thus the form of 


f required in equation (30) is 
/ men 5 OY 
The proposed strain energy function is 
- Ap _ 
which gives the following constitutive equation: 
(e)i. pine syle Y 


where y > 0. The constants which provide a reasonable fit to the 


data are 


O 
Ul 


350 psi 


= - 1420 psi 
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and this is curve D in Figure 4. Curve D follows the deformation of 
the room temperature sample very closely. To fit the cooler temper- 
ature data, a smaller value of y is required. 

In Figure 4 curve A represents the neo-Hookean prediction 
with C, = 1400 psi. Curve B is the Mooney representation using 
C, = 1400 psi and Cy = 1272 psi. Curve C is a single term version 
of Ogden's strain energy function with C, = 750 psi and = 4, 

The viscoelastic response of the material is shown in 
Figure 5. The data shown represent the most consistent results ob- 
trained in the creep tests. Using the equilibrium value of the de- 
formation from these results, appropriate values of the constants C, 
and n in equation (60) are obtained. Also, the root R is the equili- 
brium value of the deformation. Thus with R = 0.958, C, = 800 psi and 
n = 2000 1b-sec/in? equation (60) is plotted in Figure 5 for comparison. 

Since the samples are not measured at the exact test 
temperature a correction for thermal expansion may be required. Figure 
6 shows the effect of temperature on the length of the sample. The 
maximum change in length is 0.01 inch for a temperature differential 
of 70°F. Thus the coefficient of thermal expansion for this material 
is 0.000143 in/°F. In removing the sample from the environmental 
chamber, the sample temperature rises as much as 10°F before the 
measurement is made. Thus the difference between the measured length 
of a sample and its actual length at the test temperature is less 
than 0.001 inch. This difference is assumed to be negligible since 
the resolution of the micrometer is only 0.001 inch. 

Figure 7 shows the effect of the sample size on the response. 


The shorter samples appear to be slightly stiffer than the longer 
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samples. 

The use of a subpress may require a correction for the sub- 
press deformation. For the range of loads used the deformation of 
the subpress is elastic and is shown in Figure 8. The maximum de- 
formation is 0.003 inch which is insignificant when compared to a 


sample deformation of more than 0.5 inch. 


Discussion of Results 

The static elastic response of the polyurethane material 
shown in Figure 4 appears to differ from the tension results ob- 
tained with other rubber-like materials by previous investigators. 
Earlier results show a continuously increasing slope to the stress- 
strain curve as i increases. The present investigation of a poly- 
urethane material in compression indicates that an inflection point 
exists in the stress-strain curve. As A decreases from 1.0, the 
Slope of the stress-strain curve decreases until the inflection 
point is reached. Experiments performed in equi-biaxial tension, 
which is equivalent to simple compression, by previous investigators 
also show this slope characteristic and inflection point. 

While previous investigators had some success in using 
simple models to predict the elastic response, this is not the case 
in this investigation. All three of the simple models attempted are 
inadequate in predicting the response characteristics of the poly- 
urethane. The models can predict the response at small deformations 
(, > 0.95) and very large deformation (A < 0.5), but do not show an 
inflection point or the transition zone of moderate strains. The 


proposed model does predict the response curve but it lacks the 
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generality of the simpler models. It is valid for uniaxial compression 
only and the strain energy function used in its derivation does not 
meet symmetry restrictions. 

Comparison of curve D with the experimental data reveals a 
very good correlation between the proposed model and the room tempera- 
ture data. The response data obtained at cooler temperatures show a 
stiffening effect in that these curves are flatter. One explanation 
of this behavior is that the material exhibits a change in behavior 
at cooler temperatures. This explanation appears to be reasonable 
since a gradual change in the response curve is observed with the 
temperature change. However, the change is primarily at large defor- 
mation and at small deformations there appears to be no change in the 
response curve with a change in the temperature. Also barrelling of 
the sample occurs for fairly moderate strains (A < 0.9) and the 
error introduced by this barrelling is not known. Thus it is possible 
that the difference in the response data is caused by a change in 
the performance of the lubricant at the different temperatures. 

The lubricant is a silicone release agent. Other lubricants 
tried include machine oil, grease, petroleum jelly, graphite oil, and 
calcium oleate. With the exception of calcium oleate, all these 
lubricants are squeezed out upon application of some load and the 
sample barrels. Barrelling occurs even with calcium oleate but it 
occurs at larger isigstee The problem with calcium oleate is the 
difficulty in application and its significant amount of deformation 
under load. The silicone lubricant gives the best compromise of 


lubrication qualities with convenience. 
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The results of the creep tests are inconclusive since the 


scatter in the data obtained is excessive. Except for the data pre- 
sented in Figure 5, no definite trends are observable in the creep 
data. One reason for the poor data is the inaccurate step loading 
technique required by the testing machine. Since there is no pro- 
vision for a perfect step function, such a load is applied by pre- 
setting the required step electronically and then switching on the 
hydraulics to execute this loading. The resulting load function is 
inconsistent, varying with each preset load, and inaccurate (over- 
shooting the preset load). Other reasons include the lubrication 
problem and the temperature effect. 

The viscoelastic model presented in Figure 5 is based ona 
neo-Hookean spring in parallel with a linear viscous dashpot. Poor 
correlation with data is probably due to the assumption of a neo- 
Hookean spring which proved inadequate in the static elastic case. 
The use of a more complex nonlinear spring results in a differential 
equation which requires numerical methods to solve. Since the data 
is poor more complex models are not attempted at this time. 

Because of the unknown errors involved in this investiga- 
tion a quantitative error analysis is not made here. Some of the 
sources of errors include the digital voltmeter in the control panel 
of the testing machine, the load control, the micrometer, the thermo- 
stat, and the thermocouples. The digital meter has a load resolution 
of 10 pounds and a displacement resolution of 0.001 inch. The load 
control is capable of maintaining a preset load to a tolerance of 
: 10 pounds. This results in a displacement accuracy of + 0.001 


inch. The micrometer has a resoltuion of 0.001 inch. The thermostat 
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differential is 2.5°F and the resolution of the thermocouples is + 1°F. 
Of the unknown errors, the most important is probably the effect of 
barrelling. . Another is the previous history of the samples. Since 
the samples are supplied by a commercial firm, their history is not 
known exactly. Some samples appear to have more entrapped bubbles 
than others. Despite these unknowns the results obtained for the 
static case are consistent. 

Informal recovery data indicate that over 90% of the 
original length is recovered within the first ten hours after un- 
loading. Eventually samples which were loaded to 1200 pounds recover 
99% to 100% of their original length; those loaded to 1800 pounds re- 
cover 93% of their original length. This indicates that permanent set 


occurs in this material for some load between 1200 and 1800 pounds. 


Appendix A contains the data for the elastic response curves. 
The data for the viscous response curves are contained in Appendix B. 
The recovery data are given in Appendix C. Calibration data for the 
Gilmore testing machine are in Appendix D. Finally, Appendix E con- 


tains a list of the materials and equipment used. 
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CHAPTER IV 
CONCLUSIONS 


The static elastic response and the viscoelastic response 
of a polyurethane material under compression are studied in this 
thesis. For the elastic response the simple models are shown to be 
inadequate for compression deformation. _A new model is proposed 
which correlates with the data but lacks generality. The Viscoelastic 
response data are relatively inconclusive, but they do indicate that 
a more complex version of the nonlinear Kelvin-Voigt model is required. 
This conclusion agrees with the results of the elastic response tests. 

In the range of temperatures studied, the temperature does 
not appear to have a significant effect on the elastic response of 
the material. This is especially true at small deformations. At 
large deformations a stiffening effect is observed with cooler temper- 
ature. 

Thus these compression tests indicate a behavior that differs 
from that observed by previous investigation using tension tests. To 


confirm this, it is recommended that a tension test be performed on 


this material. 
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APPENDIX D 
CALIBRATION 


Calibration of Load Cell 

The load cell of the Gilmore testing machine is calibrated 
using Morehouse Proving Rings #4005 and #4006. Proving ring #4005 
has a capacity of 500 1b. compression and #4006 has a capacity of 
4000 1b. compression. The rings are calibrated by dead weights by 
National Standards Testing Laboratory. The temperature during cali- 


bration is 70°F. 


CALIBRATION DATA - RING #4005 


Machine Instrument Temperature Corrected Force 


Reading Reading Correction Reading Applied Enron 
1b. div. div. div. 1b. TD: % 
50 92.45 0.06 G7 Tol 5202 2.62 4.98 
100 212.60 OLis Ze 2/3 108.94 G2940) 8.21 
200 456.90 0.28 457.18 232.38 3230) once 
300 660.25 0.40 660.65 330.71 35.71, dahon | 
400 785.05 0.47 785.52 S95ea7 4°73, ¥s20 


a —— —— — 
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CALIBRATION DATA - RING #4006 


Machine Instrument Temperature Corrected Force 


Reading Reading Correction Reading Applied Error 

Tb. div. div. div. 1b. Tb. % 
400 102.30 0.06 102.36 413.9 13.9 3337 
500 125.70 0.08 1252/8 509.5 pas 1,87 
600 152.00 0.10 152.10 615.7 TS 7 2.50 
1000 Zoo sO 0.16 253.86 1024.5 24.5 2209 
1200 303.60 0.19 303219 pe 1224.4 24.4 1.99 
1600 404.60 0.25 404.85 1626.8 26.8 165 
2400 610.10 0.37, 610.47 2438.9 38.9 1.60 


Calibration of Displacement Transducer 

A dial deflection gauge is used to check the displacement 
transducer. The transducer is accurate to within 0.001 inch of the 
dial gauge reading over the entire 6 inch travel of the compression 


piston. An Oldak Precision Dial Indicator, model MS38,is used. 
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APPENDIX E 
LIST OF MATERIALS AND EQUIPMENT 


Materials 
Polyurethane elastomer: Vibrathane 6006 (Uniroyal Chemical) 


Silicone release agent: GE SF-96A (General Electric) 


Equipment 
50,000 1b. test frame: model 426-50 (Gilmore Industries) 
Subpress: model 85573 (Tinius Olsen Testing Machine Co.) 
Millivolt potentiometer: model 8686 (Leeds and Northrup Co.) 
Graphic recorder: model G-14A-1 (Varian Associates) 
Switching box: model 220 (Baldwin-Lima-Hamilton Corp.) 
Refrigeration unit: special model (Associated Refrigeration 

Industry) 

Micrometer: caliper type (Moore and Wright) 
Proving rings: #4005 and #4006 (Morehouse) 
Dial indicator: model MS38 (0ldak) 
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